Chain Collapse and Counterion Condensation in Dilute Polyelectrolyte Solutions 
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A new quantitative theory for polyelectrolytes in salt free dilute solutions is developed. Depending 
on the electrostatic interaction strength, polyelectrolytes in solutions can undergo strong stretching 
(with polyelectrolyte dimension R g ~ l g ' 3 N, where Ib is the Bjerrum length and N is the number 
of the chain segments) or strong compression (with R g ~ Ig 1 ^ 2 N 1 ^ 3 ). A strong polymer collapse 
occurs as a first-order phase transition due to accompanying counterion condensation. 

PACS numbers: 36.20.-r, 61.25.Hq, 64.10.+h 



Many important synthetic and biological macro- 
molecules are polyelectrolytes and their properties dif- 
fer significantly from that of the neutral polymers [Q. 
fmportance of the former stimulates a variety of their 
analytical H and numerical H studies. However, even 
for a simplest system of a salt-free solution of linear flexi- 
ble polyelectrolytes (an electroneutral system of charged 
chains and counterions) , many fundamental properties 
are still unclear or described on a phenomenological level 
only. 

The neutral chains undergo a collapse transition as the 
solvent quality decreases fuEl , but collapse and stretching 
of polyelectrolytes follow quite different laws. This hap- 
pens due to a particular role of the counterions in these 
processes, which can not be reduced only to a simple elec- 
trostatic screening and to an increase of the persistence 
length B. Under conditions of extreme dilution, for a 
weakly-charged polyelectrolytes, the counterions occupy 
the whole volume almost uniformly, with a very low con- 
centration, owing to the entropy "forces" . Under these 
conditions, the polyelectrolytes tend to strong stretching 
caused by the strong (weakly screened) inter-segment re- 
pulsions; the chain dimensions in this case are propor- 
tional to a number of charged segments. Contrary, under 
a strong polyelectrolyte charge, when the electrostatic en- 
ergy of polymer-counterion attraction is larger than the 
corresponding loss of entropy due to the counterion lo- 
calization, an essential part of the counterions localizes 
in the close vicinity of the polyelectrolytes. This effect, 
called the counterion condensation ||, leads to an effec- 
tive polymer charge screening and can cause an essential 
decrease of the polymer size. The corresponding chain 
collapse can be even stronger than for the well-studied 
neutral chains. 

In this letter we present a new quantitative theory of 
the polyelectrolytes in salt free dilute solutions. We de- 
scribe the polyelectrolyte dimension in a wide range of 
parameters, covering the area of the counterion conden- 
sation, and analyze the nature of the condensation. 

Consider a dilute solution (electroneutral as a whole) 
of charged chains, each composed of N >■ 1 segments 



(monomers) with the bond-length a. Let each monomer 
have a unit charge e (say positive) so that N counte- 
rions of the opposite charge (say negative) are present. 
Let V = (471-/3)^ be the volume per one chain with 
the R ws being the Wigner-Seitz radius, and we assume 
that the solution is very dilute, so that the conditions 
R g -C R W s and ro R ws hold. Here R g is the gyration 
radius of the polyelectrolyte chain and rjj is the Debye 
screening length built up on the counterion subsystem. 
Thus, we consider the case when the counterion clouds of 
the different polyelectrolyte chain do not overlap so that 
we can concentrate on the one-chain problem. The Bjer- 
rum length in the system is Ib = /3e 2 /e, where e is the di- 
electric permittivity of the solvent, j3 = l/(fcsT) with T 
and k.B being the temperature and the Boltzman factor, 
respectively. We also assume that the point-counterions 
approximation may be used. 

We use the following Hamiltonian for the one-chain 
problem of interest: 



H = H„ 



ch 



l.ch 



ch 
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Here -£f n . c h is the Hamiltonian of the neutral chain, which 
accounts for all noncoulombic interactions between the 
monomers of the chain, ifei.ch accounts for the Coulom- 
bic interactions between the monomers of the chain, H c 
contains the ideal counterion part and the Coulombic 
part of the counterion-counterion interactions; finally, 
H c ~ch accounts for the Coulombic interactions between 
the counterions and the chain. It is convenient to write 
the Coulombic interactions in terms of the microscopic 

N 

densities of the counterions, p c (r) = ^ 8 (r — r,*), and 

N 

monomers, p m (r) — 5 (r — Rj), where {r;} and {Rj} 

are coordinates of the counterions and of the monomers 
respectively. ff e i. c h and i? c -ch read in this notations Q: 



0H el . ch | 



drdr'ip (r - r') p m (r) p m (r') , (2) 
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-h 



drdr'ip (r - r') p m (r) p c (r') 



(3) 
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where (e 2 /e) (f(r) — (e 2 /er) is the Coulomb poten- 
tial. For the subsequent analysis it is worth to map the 
counterion part, H c onto the Hamiltonian of the one- 
component-plasma (OCP) S. The OCP model is formu- 
lated as follows: the point charges are immersed into the 
structureless compensating background of the opposite 
charge. The background charge density is ep(r), and the 
average local density of the point charges (counterions in 
what follows) is equal to p(r) [|| . The OCP Hamiltonian 
reads: 

Hocp [p(r)] = H c + H hh \p( T )] + H hc [p(r)} , (4) 

where (3H hh = {l B /2) f drdr'ip(r - r>(r)p(r') de- 
scribes the background self-interaction, while /3-ffbc = 
— Ib J drdv'<f[r — r')p(r)p c (r') gives the energy of the 
background-counterion interaction. From Eq. (Q) it fol- 
lows that H c = H OCP - H hc - H hh . 

To address the problem of the gyration radius of the 
chain, we consider the conditional free energy of the 
system, F(R g ), as a function of the gyration radius 
R g . It is convenient to define an effective Hamilto- 
nian of the neutral chain, -ff n .ch, as exp (— /3-ffn.ch) = 
D(R g ) exp (— /3-ffn.ch); where the conditional function 
D(Rg) is equal to unity if the coordinates of the 
monomers are consistent with the condition for the gy- 
ration radius to be R g , otherwise it equals zero. Since 
the free energy and structural properties of the systems 
with Hamiltonians H n c h, and -ffocp a re known (to some 
extent) , it is reasonable to write the total Hamiltonian of 
the system as a sum of two parts: the reference part, 

H = H n . ch (Rg) + Hocp[p(r)} , (5) 
and the perturbation part, 

W' = Y / drdrV(r - r') [p m (v)p m {v') - p(r)p(r')] 

-l B J dvdv'^v - r') [p m {v) - p(r)} p c (r') , (6) 

and use then the Gibbs-Bogoljubov inequality: 

F(Rg) < F n . ch (R g ) + F OGP [p(r)} + {H') Hq . (7) 

Here F n , c \,.{Rg) is the free energy of the neutral chain 
with the gyration radius R g . -Fbcp[p( r )] is the free en- 
ergy of the inhomogeneous OCP with the background 
charge density p(r). Finally, (H') H is obtained by aver- 
aging the perturbation part, if', given by Eq. (0) over 
the reference Hamiltonian, Hq, given by Eq. (§). Min- 
imizing the right-hand side of Eq. (Q) with respect to 
p(r), which is equal to the average density of the coun- 
terions, one obtains an estimate to the conditional free 
energy, F(R g ). Minimizing then F(R g ) with respect to 
R g , one finds the gyration radius. 

The free energy of the neutral polymer is a sum of 
the elastic part, written in the Flory-type approxima- 
tion jl| as ksT^ (a 2 + a~ 2 ) and the interaction part, 



written on the level of the second virial approximation 
as k B TBN 2 / (4tt R 3 g / '3) . Here a is the chain expan- 
sion factor, a 2 = R 2 g /R 2 g- lA , with R g .id being the mean- 
square gyration radius of the ideal chain, R 2 id = Na 2 /6, 
7 ~ 9/4 for the Gaussian polymers and B is the 
second virial coefficient. Here we consider the case of 
a good solvent, B > 0. With the reduced coefficient 
B* = &l 2 Bj (37ra 3 ), we write for F n . ch (R g ): 

/3F n . ch (R g )^j(a 2 +a- 2 +B*N 1 / 2 a- 3 ) . (8) 

The OCP-part of the free energy reads: 

(3F OCP [p(r)} = J drp(r)ln [A 3 c p(r) - l] 

+ J drp(r)*ggpWr)] , (9) 

where the first term in the right-hand side of Eq. (|^) 
refers to the ideal part of the OCP-free energy (A c is the 
thermal wavelength of the counterions), while the sec- 
ond term is the excess part, written in the local den- 
sity approximation (|]. The function ^ocpK 1 ")] de- 
notes the excess free energy per ion in the inhomoge- 
neous OCP-model ||. It may be expressed in terms 
of the (local) "plasma" parameter T = Ib/o-cj where 
a c = {3/[47rp(r)]} 1 / 3 is the (local) ion-sphere radius of 
the counterion. For r « 1 ( the Debye-Huckel limit) 
v|/ OC p = -V3r 3/2 , while for T > 1 an analytical fit for 
the Monte Carlo (MC) data is available j9) . One can also 
use an approximate equation of state (EOS) for the OCP, 
which covers all the range of T |To) . It gives a reasonable 
10% accuracy in the range of < T < 100, but loses con- 
siderably its accuracy for T > 100. Recently we proposed 
much more accurate approximate EOS for the OCP [ pr| : 

*ocp(F) = 7 [ln(l + cT) - cT] - ^(cT) 3 / 2 arctan -L= , 

4 i vcr 

(10) 

where c = (2/3) (2/tt 2 ) 1/3 . This EOS has the correct 
Debye-Huckel limit. It agrees within 1 — 2.5% with the 
MC data for the most of the range of T, and has a max- 
imal deviation from the MC data (about 10% ) in the 
interval 0.1 < T < 0.5. So we use Eq.© here. 

Finally, taking into account that (p c (r)) Ho — p(r), we 
write for the perturbation part: 

0{H') Ho = l fj dvdv' [<p (r - r') g 2 (r, r') (11) 
-2ip (r - r') 5 i (r) p (r') +<p(r-r')p (r) p (r')] , 

where g 2 (r,r') = (p m (r) p m (r')) Ho is the pair correla- 
tion function for the monomer- monomer density inside 
the macroion's core (i.e. in the volume confined by the 
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radius of gyration) and g\ (r) = (p m {*)) H is the average 
monomer density inside the core. 

For the counterion density distribution we adopt here 
a simplified model. We introduce two characteristic den- 
sities: pi n , the average counterion density inside the 
macroion's core, and p out , that for the outer region. This 
approximation implies that the size of the transient re- 
gion where the density of the counterions changes from its 
in-core value to the bulk value is small compared to the 
radius of gyration. From the normalization condition, 
(4tt/3) (R 3 WS - R 3 g ) pout + (47r/3)i?3p in = N, we deduce 
that pout is determined by p m . It is convenient to use 
a dimensionless density, p = pi n /n, where n = N/V g 
is the average monomer density inside the core and 
V g = (4:Tr/3)Rg is the gyration volume. The monomer 
correlation functions are approximated as g\ ~ n, and 
g2 ~ n ■ n = n . 

Using these approximations one can find all contribu- 
tions to the total free energy. For example, the pertur- 
bation part reads: 



N 



5 \R, 



2Rg 

3xt<Mj 



(12) 



where the terms 0(R 3 g /R 3 us ) are omitted. The OCP- 
part is also easily computed; it is somewhat cumbersome 
to be written explicitly in a general case. We analyze be- 
havior of the system in the limit R ws 3> R g and N 3> 1. 
Keeping for the free energy only leading terms with re- 
spect to vanishing R g /R ws and 1/iV, analyzing relative 
contribution of different addendums and omitting less 
significant ones, we finally arrive at the following result 
for the total (conditional) free energy: 



PF(a,p) 9 
N AN 



2 1 

a 2 + 

cr 



B*N 1 / 2 



or 



1/3 ^4/3 



N 1 / 6 



a 



O 1* ATl/2 

(l- P f 

5 a 



-3(l-p)lnR* ws 
(13) 



B = Z_B6 1 / 2 /a and 



where we introduced dimensionless l* E 

R-ws — Rws/o,. 

The equilibrium state of the system is determined by 
the free energy minimum with respect to both variables a 
and p. Eq. (|k|) clearly demonstrates a competition for 
the equilibrium "in-core" counterion density p between 
the two largest (at p < 1) terms, second and fourth in 
the right-hand side of Eq. (13). The second, negative 
term, large for R ws 3> 1, accounts for the counterions 
entropy It tends to minimize the free energy by mini- 
mizing p; i.e. it drives the counterions apart from the 
polyelectrolyte in order to fill uniformly all the space of 
the Wigner-Seitz cell. The positive fourth term (also 
large at N ^> 1) accounts for the free energy of screened 
Coulombic interaction between monomers. It is minimal 
if all the counterions are condensed on the polymer, i.e. 



when p = 1 • Thus mainly this two-term competition de- 
termines the equilibrium counterion density. The third 
term in the right-hand side of Eq. ([l3|) , which accounts 
for the counterion density fluctuations, becomes impor- 
tant only if p — > 1. Next, the free energy is to be min- 
imized with respect to the expansion factor a. At this 
step, the first term in the right-hand side of Eq. ( |l3| ) 
(which does not depend on p) becomes important. The 
following asymptotic cases give explicit solutions: 

(i) If l^N 1 / 2 <C a\nR^ s , the second negative 
counterion-distribution entropic term is the most impor- 
tant in Eq. Q and p -> as - {R g /R ws f. This 
case corresponds to the unscreened Coulombic interac- 
tions between the chain segments, so that polyelectrolyte 
tends to expand, and a> 1. We can write in the leading 
terms now 



(3F 
N 



9a 2 
1~N 



3 l* B N^ 2 



const 



(14) 



Thus the equilibrium expansion factor and the radius of 
gyration are, respectively, 
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1/3 



^1/3^1/2 



Ra 



l* B 1/3 N . (15) 



This regime corresponds to the polyelectrolyte strong 
stretching. 

(ii) If IgN 1 / 2 > alni2* s , the fourth positive, inter- 
segment screened- Coulombic interaction term in Eq. dl3| ) 
is essentially larger than the counterion entropic term, 
p ~ 1 and q« 1. Thus, 



(3F 



9B* 



4iVl/2 a 3 



1/3 



NV*t 



(16) 



and the equilibrium dimensions are 



3 G 



1/3 b* 1 / 2 



Ra 



B*i/2jyi/3 



* 1/2 



ATI/6 



or 



a/2 



(17) 



This is a regime of the polyelectrolyte strong collapse. 

For the general case we solved the minimization prob- 
lem numerically (again for the conditions, R ws R g 
and N 3> 1). We analyzed the dependences of R g on 
Ib at fixed B* , R^ s and N. The i?£,(^)-dependence for 
some particular values of the parameters are shown in 
Figure |l|. For small and large values of the reduced Bjer- 
rum length, l* B , the radius R q is changing in accordance 
with the asymptotic Eqs. ( [l5| ) and (|l7|). Our findings 
are in qualitative agreement with the results of the nu- 
merical study in Ref. ||, where the same two different 
regimes in the i? g (/^)-dependence were observed. The 
most interesting is however the case of the intermedi- 
ate values of l* B . In this region we observed a sharp 
bend in the dependence of the equilibrium free energy 



3 



on l* B with a discontinuity of its first-order derivatives. 
We interpret this as a first-order phase transition from 
the strong stretching regime, with a ^> 1, to the strong 
collapse regime, with a < 1. This phase transition is 
accompanied (or driven) by the process of the counte- 
rion condensation, when the counterion density changes 
from p -C 1 (counterions are uniformly spread over the 
bulk) to p ~ 1 (practically all counterions are confined 
inside the polyelectrolyte globule). After the counterion 
condensation the polyelectrolyte dimensions become es- 
sentially smaller than they would be for a neutral chain 
with the same volume interactions. 




CL 0.5 



FIG. 1. (a) The logarithmic dependences of the poly- 
electrolyte expansion factor a versus the reduced Bjer- 
rum length, l* B = IbV^/o, calculated for N = 10 3 , 
Rws — 10 5 with B* = 1. A dramatic jump-like decrease 
of the polyelectrolyte size occurs along with the first order 
phase transition due to the counterion condensation on 
the polyelectrolytes. (b) The corresponding dependences 
of the near-polyelectrolyte counterion density. Dotted 
curves refer to the numerical analysis of the simplified 
Eq.(|l3|). Solid curves account for the next-order correc- 
tions with respect to the small parameters (R g /R wa ) and 
(1/N). 



The qualitative explanation to this phenomenon fol- 
lows from the fact that for the strong Coulombic interac- 
tions (l B » 1) the OCP-part of the free energy gives rise 
to a negative pressure which takes over the ideal chain 
and counterions cntropic pressure. It may be balanced 
by the intersegment noncoulombic repulsion but only at 
some degree of compression. 

In conclusion, we developed a simple theory of the lin- 
ear polyelectrolyte salt-free dilute solutions. We analyzed 
the dependence of the gyration radius of the chain on 
the Bjerrum length, Ib, which characterizes the strength 
of the Coulombic interactions in the system and found 
two different regimes in this dependence, the same as 
were observed in numerical studies. Additionally, we de- 
tected a first-order phase transition from chain stretching 
to strong collapse, which is accompanied by the counte- 
rion condensation on the polyelectrolyte. 

Helpful discussions with B.Weyerich are highly appre- 
ciated. 
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